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Abstract 

We present a theoretical expression for the spin-echo decay rate, I/T2G, in 
the quantum-critical regime of square lattice quantum antiferromagnets. Our 
results are in good agreement with recent experimental data by Imai et al. 
[Phys. Rev. Lett. 71, 1254 (1993)] for La2Cu04. 
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I. INTRODUCTION 



Part of the recent interest in high-T^ superconductivity has been devoted to the study of 
low temperature magnetic phases of undoped and weakly doped antiferromagnetic La2Cu04 
and related compounds. Neutron scattering measurements of the correlations length at 



In this region, the correlation length is determined predominantly by classical, thermal fluc- 
tuations, while quantum fluctuations only renormalize the values of input T = parameters 
- the spin stiffness ps and the spin-wave velocity c. The RC value of the correlation length 
was also successfully used |^ to fit the data on longitudinal spin-lattice relaxation rate, 1/Ti, 



More recently, significant attention has also been devoted to a region of intermediate 
temperatures where the temperature becomes larger than the energy scale associated with 
the spin- stiffness at T = 0. Under these circumstances, the system behaves almost as if 
Ps = 0, i.e. as if it is at the quantum transition point separating the Neel and quantum- 
disordered states. This region, where the temperature is the largest energy scale in the 
problem, was first identified in [§] as the quantum-critical (QC) region, and then studied in 
detail in Clearly, in the QC region, all Bose factors are of order of unity, and hence 

quantum and classical fluctuations are equally important. 

It seems appropriate at this point to review some of the general ideas on the RC/QC 
crossover, and to present the aims of the theoretical comparisons with experiments. All 
two-dimensional magnets with a Neel ordered ground state are characterized by three energy 
scales: ps, T, and a near- neighbor exchange constant J. It was pointed out in Ref. that 
when Ps/J and T/J are small, the properties of such antiferromagnets are characterized by 
universal functions of T/p^. The physics is a smooth function of T/ps, especially simple 



T < 500K have established that it agrees very well with the theoretical result for the 
low temperature renormalized-classical (RC) region |^-^ (we set /cs = 1) 




+ ... 



(1) 



for AOOK <T < QOOK |[. 
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in two limits: the RC limit {T/ps — > 0) and the QC limit {T/ps oo). While only 
approximate N = oo results have been obtained for the full functional dependence on T/ps, 
rather precise numerical predictions are available in the two limits. At intermediate values 
of T/ps, deciding which of the two limits is appropriate is subject to interpretation, and 
different measurements may lead to different conclusions: we shall present a case for our 
choices below. A second, important issue is that of the corrections of order Ps/J and T/J. 
These corrections are in fact non-universal, depend upon the nature of the lattice cutoff, 
and impossible to calculate precisely in an analytical theory. As theorists, the best we 
can do is to calculate the universal parts, compare with the experiments and numerical 
simulations, and then determine whether the discrepancies can be attributed to these non- 
universal corrections. Clearly, Ps/J corrections will eventually become important at high 
enough temperatures - this will then impose a model-dependent upper boundary of the QC 
region. It turns out however, that in the temperature range of experimental observations, 
the universal terms agree surprisingly well with the data indicating that these non-universal 
corrections are in fact quite small. Why this should be so is not completely understood, 
but we can offer the following plausibility arguments: (z) For small T/ps, Hasenfratz and 
Niedermayer |^ have shown that there are no Ps/J corrections to the leading terms in 
the low temperature expansions of observables; (ii) numerical studies of the dynamical 
susceptibility, x{k,uj), have shown that in the whole temperature range of experiments 
it remains strongly peaked at (vr,7r), the local susceptibilities related to NMR and NQR 
measurements are then given by momentum integrals confined to A; ^ A where A is the 
upper cutoff in momentum space; {in) recent numerical analysis of the cutoff-dependent, 
two-loop equations of the sigma model [|10[, shows almost no cutoff dependence in the QC 



results until T ~ 0.5cA. We will argue below that the upper boundary of the universal 
QC region in the spin-1/2 square lattice Heisenberg model (above which T/J terms cannot 
be neglected) is about 0.6 J (~ QOOi^' in La2CuOi, we use J = ISOOi^). This temperature 
approximately coincides with the maximum temperature of experimental measurements. 
We begin by presenting our scenario of where the RC/QC crossover occurs in the universal 



functions of T / ps- Despite the fact that in Eqn. (0) the dimensionless ratio is T/2'Kps, two of 
us have recently shown exphcitly that the parameter which governs the crossover behavior 
between RC and QC regimes is three times larger: x = 3T/27rps. The reason for the extra 
factor of three is the following. From renormalization-group studies, it is known that the 
effective spin-stiffness which has to be compared with the temperature, is p^-^-^ = 2ti ps/N 
where is the effective number of components of the order parameter which contribute to 
the coupling constant renormalization Deep in the RC region, longitudinal fluctuations 
are suppressed, and the running coupling constant diverges at the scale of correlation length 
only due to interactions between — 1 transverse spin- wave modes (for 0{N) systems). 
In this situation, N = N — 2, i.e, = 1 for 0(3) case. On the contrary, in the crossover 
region near the quantum phase transition, all fluctuation modes contribute equally to the 
correlation length and we have simply N = N. For 0(3) magnets then p^jf = 2ttPs/3, and 
therefore x = T/p^ff = 3T/2ttPs. It is the variation of p^-^^ vs T/ps which eventually gives 
rise to — 2 substitution for pe// deep in the RC region, as was observed in g]. 

For the S" = 1/2 antiferromagnet on a square lattice, both perturbative [|n[] and numer- 
ical studies yield ps = 0.18J, c = 1.67Ja, so that x = T/0.38J. On general grounds, 



we expect the crossover between RC and QC regimes to occur around x ~ 1, although 
not necessarily at the same x for all observables. The upper boundary of the universal QC 
region is, we said, around 0.6J. Then, in the undoped antiferromagnets, the temperature 
range of the QC behavior should be 0.4J < T < 0.6J, which for La2CuOi corresponds to 
GOOK < T < 900K - a range which is accessible to experimental studies (a wider QC region 
is expected in the doped cuprates [0). The uniform susceptibihty data in this range IT^JTBI] 
was compared to both RC [|,|| and QC p[ formulas, and very good agreement with the QC 
result was found. Furthermore, the measured spin-lattice relaxation rate, 1/Ti, was found to 
be nearly temperature-independent above TOOK, which is consistent with the QC behavior; 
the measured constant value of 1/Ti ~ 2.7 x lO^sec"^ is also quantitatively reproduced by 
the QC theory of Ref § which predicts 1/Ti ~ (3.2 ± 0.5) x 10^sec-\ 

Very recently, Imai et al. |]16| reported results on the Gaussian component of the spin- 



echo decay rate I/T2G in the temperature range between 450i^ and 900X. Much can be 
learned from their resuhs, but we will need theoretical predictions for T2G in the QC region, 
which are computed in the following section. 



II. I/T2G IN THE QUANTUM-CRITICAL REGION 

The Gaussian component of the spin-echo decay rate I/T2G is related to the static sus- 
ceptibility, x(q) = x(q,^^ = 0), by |7 



1 \ 2 1 r d'^n 



T2gJ a' J 47r2 



d'^q ,2 



^l(q) X(q) 



(2) 



where a is the interatomic spacing and v4x(q) is a formfactor in the direction perpendicular 
to Cu02 plane, which early studies [0 estimated as A^(7r,7r) = 4.85 • WK I sec. In the 



temperature region of interest, x(q) is strongly peaked at q = (7r,7r). The second term 
in is nonuniversal in two dimensions, but it contains one less power of the correlation 
length compared to the first term, and thus can safely be neglected at low temperatures 
where ^ ^ a. To the same accuracy, we have to take the form-factor exactly at q = (vr, tt), 
in which case (1/T2g)^ measures simply the local static x^. 

Deep in the quantum-critical region, 1/Ti oc T'', oc q'"^^^ f{q/T) 0, where f{oo) 
is a constant, and rj = 0.028 |[T9|j20|l is the critical exponent for the spin correlations at 
criticality. A straightforward analysis then yields 1/^20 T~^+'''. One may combine 
this with previous results on 1/Ti, and the complete scaling forms (including prefactors) 
obeyed by the staggered susceptibility p to obtain 

TT\ _ ( A^in^y he 
T2G~[A^^in,n)) a 

where 7?. is a universal number, computable in the expansion. Here A\\{7i,it) is the 
in-plane hyperfine formfactor which appears in the 1/Ti measurements of Ref. H. The ratio 



TT1/T2G is indeed found to be temperature independent in the data of Imai et al. [116 



However for quantitative comparisons of I/T2G data with the theory, and for the estimate 
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of the correlation length, Imai et al. used a RC expression for modified by finite- 

temperature corrections in a manner first discussed by Shenker and Tobochnik [^. Below 
we take an alternative approach and compare the experimental data to the QC formula 
for I/T2G which we derive here. We will see that the agreement between our theory and 
experiment is rather good, and thus confirm the original conclusion of that the data 

on I/T2G favor QC behavior at intermediate temperatures. 

As input for our calculations, we need the expression for the static susceptibility near the 
antiferromagnetic wave vector. In the QC region, the only energy scale is the temperature, 
and the scaling function for the susceptibility depends only on g = hcq/T. Using the results 
of Ref P) , we then obtain 

where is a universal function given by 

Here Z is a rescaling factor, m is proportional (but not exactly equal) to inverse correlation 
length, and is a self-energy given by 

^'') = ^?i4^ nftl5j ■ 

where iDn = ^^n/T, Gik, ujn) = l/(k'^ + uj^ + m^), and Ii{k, oJn) is a polarization operator. 

Below we calculate I/T2G in the two leading orders in expansion for 0{N) magnets. 
The physical case will be considered at the end by setting = 3. To first order in we 
have § 

An 

A 0.23r 



Z = 1 + r/log 



2 ^2 1 ^ U.231\ 

= 02 (^i + ^log- + ^^j 



S(g)=g2(r/log^)+S(g). (7) 
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Here O = 2 log[(y5 + l)/2], A is an upper cutoff and oc 1/A^ stands for the regular 
part of tlie self-energy term. Substituting Eqn. (0) into Eqn. (^), we observe that all A- 
dependent terms disappear as they should, so that the scaling function for is universal. 



Performing then the momentum integration in Eqn. ([ 
contribution from we obtain after some algebra 

1 Al ( NTY he 



and numerically evaluating the 



2G 



av47r0 Ps V^^P 

2np - \ ^l""' 



T 



1 + 



0.22 



A^ 



NT 



(8) 



where z/ ~ 0.7 is the critical exponent for correlation length. It is then convenient to reexpress 
the result for I/T2G in terms of the actual correlation length defined from the exponential 
(e"'"/'^) decay of the spin-spin correlation function, or, equivalently, from the pole of the 
static structure factor on the imaginary q axis. From the large A^ theory of Ref P], deep in 
the QC region we have 

/27rpA^/- 



1 + 



V NT 



(9) 



Using Eqn. 



we can rewrite Eqn. 
1 Al 



as 



A^o^ f NTV ^ 



T 



2G 



1 27rp, 



0.46 



N NT 



(10) 



The advantage of using Eqn. (|T0D is in the form of the correction term which now has an 
extra factor of 1/A^; this is because at A^ = 00 all corrections related to a deviation from 
pure criticality are already absorbed into the correlation length. We will assume that the 
remaining corrections are small and neglect them below. 

We now use the values of Nq, ps and c for S = 1/2 Heisenberg antiferro magnet, the same 



value of the form-factor as in |18,21|, and rewrite the result for A^ = 3 as 



— ^ 0.546 - X 10^ sec^^. 
T2G a 



(11) 
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Finally for the ratio T1T/T2G deep in the QC region, we obtain using previous theoretical 
result H for 1/Ti 

T T 

^ 3.36 X lO^K (12) 

T2G 

Note however, that l/Ti is itself of the order and corrections to l/Ti have not 

been calculated. 

III. DISCUSSION 

We now turn to a comparison with the experimental results of Imai et al. |jT6|. We first 
use Eqn. (|ll]) and infer from the data the values of the correlation length. Fig.|l] presents our 



theoretical result obtained with no adjustable parameters together with neutron-scattering 
data available at T < SGOiT [|1| and the data of numerical simulations [0,0. We see that 
above TOOK (where 1/Ti levels off to a constant value expected in the QC regime), the 
data inferred from the T2G measurements using the QC formula are in good agreement with 
numerical results. At lower temperatures, the actual correlation length increases faster than 
in Eqn. ([Tl|) , which is a clear signature of a crossover into the RC regime. 

Next, we compare the experimental data directly with our Eqn. (^ for 1/T2g{T). The 
comparison requires some caution because the (9(0.38J/T) corrections in Eqn. (H) are not 
expected to be small. However, previous studies of 1/Ti have shown [|,^] that between 
7007^ and 9007^, the experimental data are surprisingly well described by the pure QC 
formula with no temperature-dependent corrections. We calculated I/T2G in the same way, 
and again found surprisingly good agreement with the experimental data above 700if (see 
Fig.0). At present we have no explanation why Ps/T corrections to Ti and T2G are small 
at T ~ 0.5 J; it is however difficult to be more quantitative without calculating correc- 
tions to subleading term, which has not been done. Note also that we deliberately do not 
compare the slope of 1/T2g{T) with our QC formula. The reason for not doing so comes 
from the results of numerical studies of I/T2G in a S* = 1/2 Heisenberg antiferromagnet [H. 
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The numerical results are presented in Fig.^ The two curves in Fig. ^ are the slopes of 
T2G vs the temperature, obtained using the full Eqn. (g) (these data are consistent with 
experiment for all temperatures measured), and its truncated version without nonuniversal 
corrections due to the second term and to the momentum dependence of the form-factor: 
T2 = aA~'^{J (Pq x^(g)/47r^)^^/^. We recall that these nonuniversal corrections can be ne- 
glected only if the correlation length well exceeds the interatomic spacing. We see that in 
the region of spin-echo decay measurements, T2G and T2 are quite close to each other, so 
that the absolute value of T2G (and hence ^) inferred from the experimental data should 
be consistent with the long- wavelength description. At the same time, the slopes of T2g{T) 
and T2(T) differ already by the factor of 1.7 at 900/^ which means that it is dangerous to 
compare the experimentally measured slope of T2G with the theoretical formula. 

Finally, our theoretical result for T1T/T2G, Eqn. (^2|), also agrees satisfactorily with the 
experimental value T1T/T2G ~ 4.3 x lO^i^. The difference is chiefly due to the theoretical 
result for 1/Ti, which is a bit larger than in the experiments Note also that the exper- 
imental values of T1T/T2G remain temperature independent even at smaller temperatures, 
where the correlation length already fits the RC formula. At the same time, deep in the RC 
region, one has T1T/T2G oc T^/^ PPT|]. This disagreement is not surprising as the RC result 
for T1T/T2G assumes a temperature independent uniform susceptibility Xu Hi- Numerical 



studies however indicate [T^ that Xu does not saturate until very low T < 0.2J ~ 300K. It is 
therefore likely that in the temperature range of experimental comparisons here, the highly 
nontrivial downturn renormalization of the spin-wave velocity, which leads to c(T) ~ ^/T at 
T — i> 0, does not occur and one has T1T/T2G = const even when ^ is given by Eq.(|l|). 

We now address the issue of whether it is possible to extend the QC behavior above 0.6 J. 
This issue is probably irrelevant for experimental studies, as no experiments have been done 
above 9007^^, but it is nevertheless important for the interpretation of numerical data. In 
particular, it has been recently shown that the numerical data for the correlation length at 
0.6J < T < J can be fitted by either by a modified RC expression P^ , p5 



or an mverse 



linear dependence as in the QC regime However, as has already been noted ||,|], the 
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slope of the inverse linear dependence in the QC fit was nearly twice as large as in Eqn. 
(y). Our point of view is that above T ~ 0.6 J, a mean-field description, similar in spirit 
to the N = oo approach is possible, but as we discussed in the Introduction, it should 
definitely include nonuniversal T/J terms. There are at least two sets of data which support 
the above conjecture. The first set are the data for the uniform susceptibility |T^,|TB|J2^ 



which clearly indicate that above T ~ 0.6J, the susceptibility tends to approach a broad 
maximum produced by short-wavelength fiuctuations. The second set of data, shown in 
Fig.||, are the numerical results for T2G [U- We see that above 0.6 J, not only the slopes 
but also the absolute values of T2G, T2 begin to differ substantially and their ratio reaches a 
value of ~ 1.7 at T = J. A similar nonuniversal behavior is likely to hold for the correlation 
length above 0.6J, though we cannot also exclude the possibility that the RC formula 
for the correlation length extends to higher temperatures than for other observables. The 
latter is however unlikely in view of present results and recent numerical results for doped 
antiferromagnets |p. 

To conclude, in this paper we have presented the theoretical expression for the spin- 
echo decay rate in the QC region of 2D antiferromagnets. We compared our QC result 
with the experimental data of Imai et al. |]TB|] and found good quantitative agreement in 



the temperature range between 700K < T < 900^. The temperature dependence of the 
correlation length inferred from the T2G data is in good agreement with neutron-scattering 
and numerical data. We have also argued that above T ~ 0.6J, lattice effects are relevant 
and the use of the universal low-temperature expressions for observables is unlikely to be 
justified. 
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FIGURES 

FIG. 1. The correlation length versus the temperature. Solid circles are the data inferred 



from the experiments on I/T2G' [16| using the QC formula, Eq. (|ll|). Diamonds are the neu- 
tron-scattering data of Ref. [Q]. The line represents the results of numerical studies |]l^j2^j9| at 
J = 1500K. 

FIG. 2. Experimental and theoretical results for the spin-echo decay rate 1/T2G- The solid 



circles are the experimental data of Imai et al. |16|. The asymptotic QC result, Eq. dSj), is shown 



as solid line in the QC region, T > 600i^, where the agreement with the experiment is expected. 

FIG. 3. The series expansion results Q for the spin-echo decay rates T2G, calculated using Eq. 
@ (solid line), and T2, calculated without the nonuniversal second term in (|2|) and without the 
momentum dependence of the prefactor (dotted line). 
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